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We have carried out an ab initio study of the structural, electronic and magnetic properties of 
zigzag graphene nanoribbons on Cu(lll), Ag(lll) and Au(lll). Both, H-free and H-terminated 
nanoribbons are considered revealing that the nanoribbons invariably possess edge states when 
deposited on these surfaces. In spite of this, they do not exhibit a significant magnetization at the 
edge, with the exception of H-terminated nanoribbons on Au(lll), whose zero-temperature magnetic 
properties are comparable to those of free-standing nanoribbons. These results are explained by the 
different hybridization between the graphene 2p orbitals and those of the substrates and, for some 
models, by the sizable charge transfer between the surface and the nanoribbon. Interestingly, H-free 
nanoribbons on Au(lll) and Ag(lll) exhibit two main peaks in the local density of states around 
the Fermi energy, which originate from different states and, thus, do not indicate edge magnetism. 

PACS numbers: 



Graphene [1] with its remarkable electronic and 
transport properties [2L in particular a high room- 
temperature mobility [3], is a promising material for 
applications in information technology. While perfect 
monolayer graphene has a gapless spectrum prohibit- 
ing standard transistor applications, nanostructuring 
can induce the required band gap. Recent efforts 
focused on quasi one-dimensional graphene nanorib- 
bons (GNRs) |H3 

and zero-dimensional graphene 
quantum dots (GQDs) [8l4Io| and, indeed, revealed 
a transport gap, e.g. for GNRs with widths below 
10 nm. Theoretically, unsupported GNRs and GQDs 
with zigzag edge geometry possess spin polarized edge 
states with ferromagnetic order along the ed ge [Ill4l3l . 
Several experimental studies provide direct [10|, Il4i4l6| 
or indirect [ItJ evidence for the presence of the edge 
states also in supported GNRs, albeit without probing 
its magnetism. Such edge states might be exploited 
for a multitude of spintronics applications 0, fl8l - [2ci| , 
but so far it is unclear, if the edge states contribute to 
the measured transport properties of nanostructures at 
all. Thus, it is important to elucidate the role of the 
substrate and of edge termination. 

Recently, several groups including us investigated in-situ 
prepared GQDs with exclusive zig-zag-edges, which 
are supported by Ir(lll) [2~0-[27]. In particular, we 
revealed the absence of edge states by a combined Den- 
sity Functional Theory (DFT) and scanning tunneling 
microscopy (STM) study (28[. Our finding was explained 
by a hybridization of the graphene 7r-orbitals with an Ir 
5d z 2 surface state at Ep, which gradually decreases in 
strength from the edge towards the center of the GQD 
and, thus, prohibits a simple shift of the edge state 
towards the interior of the GQD [H|. In contrast, STM 
studies of Tao et al. [l6[ provided convincing evidence for 
the presence of edge states in GNRs chemically prepared 



from carbon nanotubes by so called unzipping [29] and 
deposited on Au(lll). A peak or a double-peak close 
to Ep was observed in the local density of states (DOS) 
close to the edge with a peak distance scaling with the 
width of the GNR. The double-peak was present for 
chiral angles up to 16.1° with respect to the zigzag 
direction and was ascribed to an antiferromagnetic 
coupling between opposite edges 16, [30[ as further 
evidenced by comparison with results from a Hubbard 
model Hamiltonian [16|. 

Here, we present a more realistic description of this 
system using DFT and including the Au(lll) surface. 
Large models were used, which enabled us to study 
GNRs of close to realistic widths and to correctly 
describe the lattice mismatch between graphene and 
Au(lll). For comparison, we also investigated GNRs 
deposited on the (111) surface of the other two group 
IB metals, Cu and Ag. For all the substrates, we 
considered both H-free and singly H-terminated GNRs. 
For H-free GNRs on Au(lll), which is most relevant 
to the experiments by Tao et al. [16[, we considered 
the two chiral angles = 0° and 5°, while only perfect 
zigzag GNRs were studied for the other systems. We 
show that all the GNRs studied exhibit edge states, but 
that the interaction with the substrate strongly depends 
on whether the GNR is H-terminated or not and, to a 
lesser extent, also on the type of substrate. As a result, 
it turns out that the edge states are magnetic only when 
H-terminated and deposited on Au(lll). Exactly this 
system was prepared experimentally recently [29], l3lj . 
Interestingly, the nonmagnetic H-free GNR on Au(lll) 
also exhibits a double peak around E^ (albeit with a 
much larger splitting as compared to experiments [l6|), 
but the two peaks originate from a hybridization of Au 
d-orbitals with the graphene and the unoccupied edge 
state, respectively, and not from a magnetic splitting of 
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TABLE I: Minimum distance between the edge C atoms of 
the GNR and the atoms of the relevant (111) surface and 
maximum distance between the center of the GNR and the 
surface. Distances are in Angstrom. 



the edge state. 

We considered GNRs with zigzag edges and a width of 
8 graphene unit cells. We employed different super cells 
in order to account for a) the lattice mismatch between 
graphene and the (111) surfaces of Cu, Ag and Au and 
b) chiral angles. Here, we describe the models having 
= 0° (the model with finite is discussed in the Sup- 
plement jlH). For Ag and Au a supercell with a size 
parallel to the GNR six times as large as the unit cell of 
the (111) surface was used. For Cu, an extremely large 
supercell would be required to account for the mismatch, 
making the calculations unfeasible. Therefore, a slightly 
compressed (about 3.8 %) Cu lattice was instead used to 
make graphene and Cu(lll) commensurate. A four-layer 
slab was used for Ag and Au surfaces, whereas the use 
of a smaller supercell for Cu (in the direction parallel to 
the GNR) enabled us to employ thicker slabs containing 
up to 12 layers. The latter calculations showed that 4- 
layer slabs are sufficient to describe semi-quantitatively 
the interaction between the GNR and the Cu surface (see 
Supplement), making us confident that the same holds 
true for Ag and Au substrates. 

The structural optimization and the calculation of the 
electronic and magnetic properties were carried out us- 
ing the plane- wave package Quantum-Espresso [33[ . Ad- 
ditional computational details are provided in the Sup- 
plement [32|. 

In the first part of the paper, we present our results 
about the GNRs on Au(lll). Ag and Au substrates are 
discussed in the second part and in the Supplement [32[ . 
Upon DFT relaxation, the H-free GNR with (9 = 0° 
on Au(lll) bends considerably: the distance between C 
atoms and the surface is much shorter at the GNR edge 
than in the interior of the ribbon (Fig. 1(a)). The max- 
imum distance between C atoms and the Au (as well as 
Cu and Ag) surface at the center of the GNR and the 
minimum distance at the edge are shown in Table I. A 
small corrugation (about 0.17 A) along the GNR edges 
in accordance with the Moire periodicity is observed. In 
particular, 4 edge C atoms out of 7 are in a quasi on-top 
configuration, whereas the remaining edge atoms are in 
a quasi bridge or hollow configuration. In the following, 
the z axis will be taken perpendicular to the surface and 
the x axis will be taken parallel to the GNR. Since the 
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FIG. 1: Structural and electronic properties of a H-free GNR 
on Au(lll). (a) Topography of the relaxed model. For the 
sake of clarity, only the top Au layer is shown. C and Au 
atoms are labeled by numbers indicating different chemical 
environments and used in (b)-(g). The top and hollow ad- 
sorption sites of the edge C atoms are denoted with the sub- 
scripts "t" and "h" respectively, (b) (p y + p z )-PDOS of a C 
atom (C m ) in the middle part of the GNR and PDOS of the 
d states of an Au atom (Au m ) located beneath, (c) PDOS of 
d states of several Au atoms starting from an atom below the 
edge of the GNR (Am) towards an atom below the centre of 
the GNR (Au 4 ). (d), (e) (p y + p 2 )-PDOS of C atoms at the 
center of the on-top (d) and hollow (e) region, in increasing 
distance from the edges; row 1 denotes the edge row. (f), (g) 
(Py + Pz)-PDOS of an edge C atom at on-top site, Cti (f), 
and at hollow site, Chi (g), and PDOS of the d states of the 
nearest neighbor Au atom (Aus and Au6 respectively), (h), 
(i) Plots of a charge isosurface of a state contributing to the 
peak at -0.2 eV below Ef (h) and at 0.3 eV above Ef (i). 



GNR is not parallel to the surface at the edge, the or- 
bitals of C forming edge states are linear combination of 
2p z and 2p y orbitals. For the same reason, the dangling- 
bond orbitals of the edge C atoms of the H-free GNR 
(which form a bonds with Au atoms upon deposition on 
the surface, see below) are also combinations of these or- 
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bitals. In Fig. 1(b) the (p y + p 2 )-PD0S of a C atom 
in the middle part of the GNR and the PDOS of the 5d 
states of an Au atom located beneath are shown: the 
contribution of the PDOS of the d states of Au near Ep 
is very small and the interaction between the two atoms 
is basically negligible. On the other hand, a comparison 
of the PDOS of the d states of several Au atoms starting 
from an Au atom below the edge of the GNR towards 
an Au atom below the centre (Fig. 1(c)) shows that, at 
the edge, the PDOS displays some peaks near Ep. Also, 
the (pj, + p 2 )-PD0S of the edge C atoms (Figs. l(d)-(e)) 
exhibit two main peaks near Ep. To shed light on the 
interaction between edge C atoms and Au atoms, it is 
useful to compare the respective PDOS at on top sites 
(Fig. 1(f)) and hollow sites (Fig. 1(g)). The (p y + p*)- 
PDOS of a C atom at an on-top site displays two peaks at 
-0.2 eV below and 0.3 eV above Ep respectively. Inspec- 
tion of the PDOS of the nearest neighbor Au 5d orbitals 
indicates that a strong hybridization between these states 
and C states occurs. In particular, a relatively high peak 
below Ey is present, in correspondence with the small 
peak of the C p states, whereas a less pronounced peak 
is found right above Ep, corresponding to the second, 
large peak of the C p PDOS. In the hollow case, two 
peaks are present in the PDOS of C p states as well, 
whereas the PDOS of the nearest neighbor Au d states 
exhibits no peak above Ep. These findings suggest that 
the two peaks in the PDOS of C p orbitals near Ep have 
a different origin. In fact, the plots of charge isosurfaces 
of two states contributing to said peaks (Figs. 1(h)- (i)) 
unequivocally show that the one above Ep is due to the 7r 
edge states of the GNR, whereas the one below Ey is due 
to states originating from a bonding between C and Au 
atoms. More precisely, the latter peak corresponds to an- 
tibonding p-d states, whereas the corresponding bonding 
states have much lower energies (about -5 eV below Ep). 
In the on-top case, the 5d z 2 and 5d yz orbitals of the Au 
mostly contribute to the bonding with C p states; in the 
hollow case, 5d zx orbitals contribute as well. Au 6s and 
6p orbitals do not play an important role in this bond- 
ing. Since the antibonding a states have lower energies 
than the edge states, the latter states are mostly unoccu- 
pied and no significant edge magnetism occurs (less than 
3 • 10 -3 /i£ per edge C atom). 

The model of the GNR with = 5° is presented in the 
Supplement. This model also bends considerably and 
exhibits essentially non-magnetic edge states for similar 
reasons (see Fig. SI). 

Next, we consider H- terminated GNRs on Au(lll). In 
this case, the interaction between the GNR and the sur- 
face is weak: the minimum distance between the edge 
C atoms and the surface is 3.45 A (see Table I). The 
GNR is almost flat, as evident from Figs. 2(a)-(b). As a 
result, the electronic properties of the GNR are weakly 
affected by the presence of the substrate: the GNR dis- 
plays a magnetic edge state with antiferromagnetic cou- 




FIG. 2: Structural and electronic properties of a H- 
terminated GNR on Au(lll). (a)-(b) Top and side view of 
the relaxed model. Only the top Au layer is shown, (c) Top 
and side view of an isovalue surface of the edge state spin den- 
sity of the deposited GNR. The system has antiferromagnetic 
order across the GNR. The red (blue) surface indicates spin 
up (down) density, (d) Non-spin-polarized PDOS of the 2p z 
orbitals of an edge C atom of the deposited GNR (Cd). The 
corresponding PDOS of the edge C atom of the free-standing 
GNR (Cf) is also shown for comparison, (e) Spin-polarized 
PDOS of the 2p z orbitals of a C atom at the left edge of the 
magnetized GNR shown in Fig. (c). 



pling between the edges and the magnetization per edge 
C atom is about 0.22 hb (comparable to the value ob- 
tained for free-standing GNRs using equivalent k-point 
meshes, see Supplement). An isovalue surface of the edge 
state spin density is shown in Fig. 2(c), whereas the spin- 
unpolarized and spin-polarized PDOS of the 2p z orbitals 
of edge C atoms are shown in Figs. 2(d) and (e) respec- 
tively. The energy splitting between spin majority and 
minority 2p z peaks on an edge is about 0.7 eV (Fig. 2(e)). 

Since the edges of the GNRs investigated experimen- 
tally in Ref. 1 16l exhibit a pronounced edge curvature, they 
resemble the H-free case of our calculation. A compar- 
ison of the experimentally observed peak splitting with 
the calculated one is tempting, although the experimen- 
tal GNR width of 8 nm to 20 nm is much larger than in 
our models (1.6 nm). The experimentally observed en- 
ergy difference ranges from 50 meV to 20 meV [l6[ and 
is an order of magnitude lower than the one from the 
DFT calculation. Thus, a direct assignment of the calcu- 
lated splitting to the experimental results is not possible, 
maybe due to a different edge chemistry, but our findings 
suggest that the interpretation of a double peak around 
Ep alone as a sign for edge magnetism might be mislead- 
ing. 

Naively, one could expect that the properties of the 
GNRs on Ag(lll) would be almost identical to those dis- 
cussed above. Indeed the electronic structure of H-free 
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FIG. 3: Structural and electronic properties of a H- 
terminated GNR on Ag(lll). (a)-(b) Top and side view of the 
relaxed model. Only the top Ag layer is shown, (c) Isovalue 
surfaces of the difference between the total charge of the GNR 
plus substrate system and the charge of the isolated (bent) 
GNR and Ag(lll). The red (blue) color indicates accumula- 
tion (depletion) of charge, (d) PDOS of the 2p z orbitals of 
an edge C atom of the deposited GNR (Cd). The non-spin- 
polarized 2p z PDOS of the edge C atom of the free-standing 
GNR (Cf) is also shown. 



GNRs on Ag(lll) is qualitatively similar to that of the 
nanoribbons deposited on Au, as discussed in the Sup- 
plement. On the other hand, H-terminated GNRs on 
Ag(lll) exhibit remarkable differences with respect to 
Au(lll), which originate from the relatively strong in- 
teraction between the GNRs and the Ag substrate. In 
fact, the minimum distance between the C atoms and 
Ag(lll) is 2.69 A at the edge (Table I) and a signif- 
icant bending of the GNR occurs (Figs. 3(a)-(b)). The 
interaction between Ag atoms and edge C atoms is partly 
covalent and a net charge transfer from the surface to the 
GNR also occurs, as evident from Fig. 3(c). In total, the 
GNR has acquired a charge of 1.6 • 10 -2 electrons per C 
atoms. This behavior is in contrast to the case of Au(lll) 
and can be qualitatively explained by the fact that Ag 
is more prone to donate electrons than Au. The PDOS 
of the C 2p z orbitals exhibits a large peak (Fig. 3(d)), 
which corresponds to the edge states. However, due to 
the n doping of the GNR, the peak is not exactly at 
but is shifted slightly downwards in energy. As a result, 
the DOS at Ep is small and no Stoner instability oc- 
curs: hence, this system is non magnetic. In fact, the 
calculated magnetization is less than 10 _3 /iB per edge C 
atom. Furthermore, the height of the PDOS peak is also 
reduced with respect to the free-standing case, due to the 
interaction with the Ag atoms (Fig. 3(d)). 

GNRs on Cu(lll) display significant structural differ- 
ences with respect to the Au and Ag substrate, which 
stem from the stronger reactivity of Cu atoms. In par- 
ticular, in the H-free case, the interaction between the 
GNR and the Cu substrate is not negligible even in the 
center of the GNR and the maximum distance between 



the latter and the surface is only 3.31 A (see Table I 
and Fig. S2a). The edge state interacts more strongly 
with the Cu substrate than with Au(lll) and Ag(lll): 
nevertheless, the PDOS bears a qualitative resemblance 
to that of the latter models (see Supplement). In the li- 
ter minated case, the strong interaction between C and Cu 
atoms also leads to shorter equilibrium distances between 
the GNR and the substrate (Table I and Figs. S5(a)-(b)). 
Charge transfer leads to the filling of the edge state in this 
model too, which is thereby non-magnetic (Figs. S5(c)- 
(d))- 

In conclusion, our simulations based on DFT indicate 
that zigzag GNRs deposited on Cu(lll), Ag(lll) and 
Au(lll) all possess edge states but do not exhibit signifi- 
cant edge magnetism, with the exception of H-terminated 
GNRs on Au(lll), whose zero-temperature magnetiza- 
tion is comparable to that of free-standing GNRs. These 
results are explained by the different interaction and 
charge transfer between the GNRs and the substrates. 
Only in the case of H-terminated GNRs on Au(lll) is 
the interaction sufficiently weak so as not to affect the 
electronic and magnetic properties of the edge states sig- 
nificantly. Hence, our simulations strongly suggest that 
experimental investigations on edge magnetism in GNRs 
deposited on metallic substrates should focus on the lat- 
ter system. 
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